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' Abstract 

o 

' The probability for quantum creation of an inflationary universe with a pair of black 

holes is computed in a modified gravitational theory. Considering a gravitational action 
, which includes a cosmological constant (A) in addition to a B? and 5 R^^ terms, the 

probabilities have been evaluated for two different kinds of spatial sections, one accom- 
modating a pair of black holes and the other without black hole. We adopt a technique 
(— I ^ prescribed by Bousso and Hawking to calculate the above creation probability in a 

semiclassical approximation with Hartle-Hawking boundary condition. Depending on 
■ the parameters in the action some new and physically interesting instanton solutions 

' are presented here which may play an important role in the creation of the early uni- 

verse. We note that the probability of creation of a universe with a pair of black holes 
is strongly suppressed with a positive cosmological constant when 5 = for a > 
but it is more probable for a < — g^. It is also found that instanton solutions are 
allowed without a cosmological constant in the theory provided 5 < Q. 
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It is now understood from the information gathered from high redshift surveys of su- 
pernovae to Wilkinson Microwave Anisotropy Probe ( WMAP ) data that the present 
universe is accelerating [1] . The universe might have also emerged from an inflationary era 
in the past predicted from COBE. In the last two decades a number of literature appeared 
which explored early universe with an inflationary phase [2]. Inflation not only opens up 
new avenues in the interface of particle physics and cosmology but also solves some of 
the outstanding problems not understood in the standard bigbang cosmology. Most suc- 
cessful theory developed so far is based on the dynamics of scalar field with a suitable 
potential which drives inflation. In order to explain present acceleration in the universe, 
various cosmological models have been proposed with exotic fields that appeared with a 
new gravitational physics considering theories other than scalar field, e.g., chaplygin gas 
[3], phantom fields [4] etc. It may be pointed out here that an alternative theory with 
a modification to the Einstein-Hilbert action has been employed to obtain early inflation 
instead of a scalar field. It is known that a modification to the Einstein-Hilbert action with 
higher order terms in curvature invariants that are become effective in the high curvature 
region lead to modifications of the standard cosmology admitting a de Sitter universe at 
early times [5]. In the same way it is important to look for modifications of the Einstein 
gravitational action with terms that might be important at extremely low curvature region 
to explain the present cosmic acceleration. Recently, various cosmological models are pro- 
posed to explain the cosmic speed up with a modification of the Einstein Hilbert action. 
Carroll et al. [6] considered a gravitational action of the form 



where R is the Ricci Scalar, = SttG, and is a mass scale of the order of Hubble scale, to 
construct such models of the universe. In the literature [7] a number of cosmological models 




(1) 
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are appeared with the above action to construct an alternative theory of dark energy and 
dark matter. However, in the theory de Sitter space is always unstable. It has been shown 
[8] that adding s, a B? term to the action allows a stable de Sitter space if the size of the 
baby universe is bigger than a size determined by ji^ in the action. Further, the theory 
containing both the 1/R and i?^ terms has better chance of passing solar system tests 
than the theory based on simple 1/R corrections [9]. It is also important to explore the 
cosmological issues in the framework of new theories which are important at the present 
epoch. The astrophysical objects, for example black holes need to be investigated in these 
theoretical framework. The mass of these objects may be greater than the solar mass or 
even less. It is known in stellar physics that a blackhole is the ultimate corpse of a collapsing 
star when its mass exceeds twice the mass of the Sun. Another kind of black holes are 
important in cosmology which might have formed due to quantum fluctuations of matter 
distribution in the early universe. These are termed as topological blackholes having mass 
many times smaller than the the solar mass. In particular, black holes whose mass (~ lO^^g') 
have life times comparable to the age of the universe so that they would be evaporating 
at present which is Hawking radiation [10]. It may provide sources of ultra high energy 
cosmic rays including extragalactic gamma-rays [11] and other astrophysical phenomena 
whose mechanisms have not been solved completely. In addition to that, the relic black hole 
might contribute to the energy density of the present universe and solve the dark matter 
problem [12]. Bousso and Hawking [13] (in short, BH) calculated the probability of the 
quantum creation of a universe with a pair of primordial black holes (in short, PBH) in (3 + 
1) dimensional universe. To compute the probablity measure, BH considered two different 
Euclidean space-time : (i) a universe with space-like sections with - topology and (ii) a 
universe with space-like section with x 5^- topology, as is obtained in the Schwarzschild- 
de Sitter solution. The first kind of spatial structure describes an inflationary (de Sitter) 
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universe without black hole while the second kind describes a Nariai universe [14], an 
inflationary universe with a pair of black holes. Considering a massive scalar field which 
provided an effective cosmological constant for a while through a slow-rolling potential 
(mass-term) BH studied the problem. Chao [15] studied the creation of a primordial 
black hole and Green and Malik [16] studied the primordial black holes production during 
reheating. Paul et al. [17] following the approach of BH studied the probability of creation 
of PBH including i?^-term in the Einstein action and found that the probability is very 
much suppressed in the i?^ -theory. Paul and Saha [18] studied probability of creation of a 
pair of black hole with higher order Lagrangian i.e., considering higher loop contributions 
into the effective action that are higher than quadratic in R. Recently, the probability for 
quantum creation of an inflationary universe with a pair of black hole in a gravitational 
theory which includes a cosmological constant (A) in addition to {6 R~^) term is evaluated 
[19]. In this paper we investigate the probability for primordial black hole pair creation in 
a modified gravity by adding an extra term aR^ to the action given by (1) in the presence 
of a cosmological constant. The above theory is important as it may be regarded as the 
synthesis of the theories one needs for early and late inflation, although it remains to be 
understood how the universe transits from early inflating phase to the present accelerating 
phase. 

Consequently it is also important to study the effects of these terms to study quantum 
creation of a universe with a pair of PBH. We calculate the probabilities for the creation 
of a universe with two types of topology namely, R x - topology and R x x S'^ - 
topology, where the later accommodates a pair of primordial black holes. To calculate 
the probabilities for these spatial topologies, we use a semiclassical approximation for the 
evaluation of the Euclidean path integrals. The condition that a classical spacetime should 
emerge, to a good approximation, at a large Lorentzian time was selected by a choice of 
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the path of the time parameter r along the axis from to 2^ and then continues along 
the T*™ axis. The wave-function of the universe in the semiclassical approximation is given 

by 

^o[hij, ^dM] « ^n^--^" (2) 
n 

where the sum is over the saddle points of the path integral, and J„ denotes the corre- 
sponding Euclidean action. The probability measure of the creation of PBH is 

P[hab, ^du] ~ 

where hab is the boundary metric and I^^ is the real part of the action corresponding 
to the dominant saddle point, i.e. the classical solution satisfying the Hartle-Hawking 
(henceforth, HH) boundary conditions [20]. It was believed that all inflationary models 
lead a density parameter to ~ 1 to a great accuracy. This view has been modified after 
it was discovered that there is a special class of infiaton effective potentials which may lead 
to a nearly homogeneous open universe with < 1 at the present epoch. Cornish et al. 
[21, 22] studied the problem of pre- inflationary homogeneity and outlined the possibility 
of creation of a small, compact, negatively curved universe. We show that a universe with 
S'^-topology may give birth to an open inflation similar to that obtained in Ref. 23. 

The paper is organised as follows : in sec. II, the gravitational action for a higher 
derivative gravity is written and gravitational instanton solutions are presented and in sec. 
Ill, the action is used to estimate the relative probability for the two types of the universes 
one with PBH and the other without PBH. Finally, in sec. IV we give a brief discussion. 

II. GRAVITATIONAL ACTION AND INSTANTON SOLUTIONS WITH 
OR WITHOUT A PAIR OF PRIMORDIAL BLACK HOLE : 

We consider Euclidean action which is given by 

lE = -r^ [ d^x^ f{R) - ^ / d\^^h Kf{R) (3) 

iOTT J OTT JdM 
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where g represents the 4-diinensional Euchdean metric, a and 5 are dimensional parameters, 
A represents the 4-dimensional cosmological constant and f{R) is a polynimial function in 
Ricci Scalar R. The second term in the action is the gravitational surface term where hij 
is the boundary metric and K = K^^Kij is the trace of Kij , the second fundamental form 
of the boundary dM in the given metric g. The second term is the contribution from r = 
back in the action. It vanishes for a universe with topology S^, but gives a non vanishing 
contribution for topology x S"^. 

(i) /S^-Topology, the de Sitter spacetime : 

In this section, we first derive the field equations from the Eucliden action and then 
explore instanton solutions in the modified gravitational theories. We now look for a 
solution with spacelike section having -topology. The corresponding four dimensional 
metric ansatz in this case is given by 



ds^ = dr^ + a^(r) dxi + sin^xi dO^ 



(4) 



where o(t) is the scale factor of a four dimensional universe and dO^ is a line element on 
the unit two sphere. The scalar curvature is given by 



R = -6 



where an overdot denotes differentiation with respect to r. We rewrite the action (3), 
including the constraint through a Lagrangian multiplier /? and obtain 



Ie 



IT 



f{R)a^-P\R + 6- + 6 



dr-^ [ d^xVhKf'(R). (5) 
ovr JdM 



Varying the action with respect to R, we determine 



P = a^f'{R). 



(6) 



Substituting the above equation, treating a and R as independent variables we get 



Ie = -^ \a^f{R) - f'iR) (a^R - Sad^ - 6a) + Qa^aRf"{R) 



dT 
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aa^f'iR) 



T=0 



(7) 



here we have ehminated d term in the action by integration by parts. The field equations 

are now obtained by varying the action with respect to a and R respectively, giving 



f"{R) 



R + 6- + 6 



d2 



1 



2f"{R)ie + 2f"{R) 



R + 2-R 
a 



+ f{R) 



0, 



4- + 2— - — + i? 



(8) 



/(i?) = 0. (9) 



We now consider f{R) = i? — 2A + aR? + ;| to solve Eqs. (8)- (9). An instanton solution 



is obtained which is given by 



o = — sin Ht 
H 



(10) 



where R = 12H'^ and H can be determined in terms of the parameters S,A from the 
constraint equation which emerged from eq.(9) : 



48H'^ - 16AH^ + 6 = 0. 



Ill) 



We note following from the Eq. (11) : (i) A real Hubble parameter is permitted for 
S < ^ , (ii) two classes of instantons are allowed for < (5 < ^ , (iii) Instanton solution 



with H 



A ; 



is obtained when S = ^4- , (iv) Instanton solution with H 



A : 



is obtained 



when 6 = 0, and (v) Instanton solution is permitted even with A = when 6 < 0. 

It is evident that the instanton solution (10) obtained here satisfies the HH no boundary 
conditions viz., a(0) = 0, d(0) = 1. One can choose a path along the r^^ axis to r = 
the solution describes half of the Euclidean de Sitter instanton S^. Analytic continuation 
of the metric (4) to Lorentzian region, xi — > ^ + ia, gives 

ds'^ = dT^ + a^(r) [-dc7^ + cosh^ a d^l 



which is a spatially inhomogeneous de Sitter like metric. However, if one sets t = it and 
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cr = + X, the metric becomes 

ds^ = -dt'^ + b^{t) [dx^ + sinh^ x ^^^^2] (12) 

where b{t) = — i a{it). The hne element (12) now describes an open inflationary universe. 
The real part of the Euclidean action corresponding to the solution calculated by following 
the complex contour of r suggested by BH is given by 

'l728aH^ + lUH'^ - 24AH'^ + 6' 



TRe _ ^ 



(13) 



With the chosen path for r, the solution describes half the de Sitter instanton with S'^ 
topology, joined to a real Lorentzian hyperboloid of topology x S^. It can be joined 
to any boundary satisfying the condition agu > . For qqm > , the wave function 
oscillates and predicts a classical space-time. We note the following cases : 

(a) solution obtained by BH [13] is recovered when 5 = 0, a = 0. 

(b) = 0, a 7^ 0, corresponds to the instanton solution obtained by Paul et al. [17]. 

(c) 5 7^ 0, a = 0, A 7^ 0, corresponds to the solution obtained by us recently [19] and 
the action reduces to / = —--^= for a vanishing cosmological constant. It is evident here 
that a realistic solution is obtained with a negative 5 only. 



3 a+ ^ 



3 5 



But 



(d) 5 7^ 0, a / and A = 0, one obtains action which is / = —4 tt 
the above solution also demands (5 < 0. 

(e) 5 = we get the action which becomes I = —12 ira — ^ . It is interesting if 
a < 0. 

(ii) Spacelike sections x S'^, the Nariai spacetime: 

We now look for a solution with spacelike sections x S"^ as this topology accommodates 
a pair of black holes. The corresponding ansatz for (1 + 1 + 2) dimensions is given by 

ds^ = dr^ + (r) dx^ + 6^ (r) dn^ (14) 
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where a(r) is the scale factor of S'^-surface and 6(r) is the scale factor of the two sphere. 
The metric for the two-sphere is given by the metric 



The scalar curvature is given by 



R = - 



d b 1\ hb 



(15) 



The Euclidean action (3) becomes 



Ie = - 



IT 



b (ii 



f{R)ab'' - /? i? + 2- + 4- + 4— + 2 



b }? 2 



b ab b"^ 62 



dT 



-— [ Vhd^xKf{R). 

StT JdM 



(16) 



One can determine (3 as is done before and obtains 



= -- f 
2 Jr 



TdM 
T=0 



/(ij)-m(H-45^-2|-|)+2mii(^+2[ 



ab^dr - TT \{ab^ + 2abb) f'{R) 



T=0 



(17) 



Variation of the action with respect to R, a and b respectively are given by the following 
equations 



2f"\R)R^ + 2f"{R) 



a b ab b^ 



R + 2r\ 
b 



0, (18) 

f{R) = 0, (19) 



2f"{R)R^ + 2f"{R) 



■ a b\ •• 



+ f{R) 
- f{R) = 0. 



„ a b ab 
R + 2- + 2- + 2— 
a b ab 



(20) 



Let us now consider f{R) = R—2A + a R^ + -^. Eqs. (18)-(20) admit an instanton solution 



which is given by 



a = — sin{HoT), b = H^^, 



H. 



R = AHl (21) 
where Hq satisfies the constraint equation 

im^ - IQKHl + 3(5 = 0. (22) 

We note the following : (i) when < 5 < two classes of instanton solutions are 



permitted, (ii) (5 = 0, instanton solution with = A is obtained, (iii) 5 = leads to 



an instanton solution with = ^, (iv) A = 0, permits instanton with = y — % when 
S<0. 

The above instanton solution (21) satisfies the HH boundary conditions a(0) = 0, d(0) = 
1,6(0) = bo, b{0) = 0. Analytic continuation of the metric (14) to Lorentzian region, i.e., 
T ^ it and x — ^ | + w yields 

ds^ = -dt^ + c^{t) da^ + dnj, (23) 

where c{t) = — i a{it). In this case the analytic continuation of time and space do not give 
an open inflationary universe. The corresponding Lorentzian solution is given by 

a(r-^"')LKe = H-^ cosh HoT^"", 

6(r^-)|,.e=^=i7-i (24) 

Its space like sections can be visualised as three spheres of radius with a hole of radius 
b = punched through the north and south poles. The physical interpretation of the 
solution is that of two - spheres containing two black holes at opposite ends. The black 
holes have the radius which accelerates away from each other with the expansion of 
the universe. This describes half of a Lorentzian Nariai universe. The real part of the 
action can be determined now following the contour suggested by BH [13], and it is given 

by 



TRe _ ^ 
-'51x52 - 



6AaH^ + 16H^ - 8AH^ + s] . (25) 
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The solution given by eq.(24) describes a universe with two black holes at the poles of a 
two sphere. It may be pointed out here that the contribution of the integrand in the action 
(17) for the instanton vanishes and the non-zero contribution of the action here arises just 
from the boundary term only. We note the following : 

(a) 5 = 0, a = 0, corresponds to the action obtained by BH [13]. 

(b) 5 = 0, a 7^ corresponds to / = —Sira — J, it reduces to a case discussed by Paul 
et al. [17]. 

(c) (5 7^ 0, a = 0, A = 0, the action is / = — ^-i==, it is real for 5 < , which 
corresponds to the result obtained by us [19]. 

(d) (5 7^ 0, a 7^ and A = corresponds to an action which is 

rRe _ Stt / 2 \ 2 

In this physically realistic solution is obtained for S <0. 

(e) 6 = the action becomes 

rRe _ ^vr/ 1\_2 ., 

III. THE PROBABILITY FOR THE PRIMORDIAL BLACK HOLES PAIR 
CREATION : 

In the previous section we have calculated the actions for inflationary universe with or 
without a pair of black holes. We now compare the probability measure in the two cases. 
The probability for creation of a de Sitter universe is determined from the action (7). The 
probability for nucleation of an inflationary universe without a pair of Black holes is given 
by 



Pgs ~ eH J. (26) 

However for an inflationary universe with a pair of black holes the corresponding probability 
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of nucleation can be obtained from the action (17). The corresponding probability is 



24aH^+4H'^-A 
^ 



(27) 



For simplicity to compare the probabilities we consider spacial cases : 

(a) The probabilities obtained by Bousso and Hawking [13] is recovered when both 6 
and a are set zero. The probabilities are 



P53 ~ e A , Pgiy^g2 ~ e A . 



(28) 



In this case the probability for a universe with PBH is less than that without PBH. The 
solution cannot admit a negative cosmological constant. 

(b) The probabilities evaluated by Paul et al. [17] are recovered here with 6 = 0. The 
probabilities are 

(29) 



-Tgs ~ e A ' J jr^giy^g2 ~ e A ' 



In this case with a positive cosmological constant and a > 0, the de Sitter universe is 
more probable. However for a < — the probability for a universe with PBH is more. 
However, the case a < leads to a classical instability in R'^ theory. 

(c) For 6^0 and a = 0, the solutions obtained in Ref. 19 are recovered. There 
are two branches of solutions admitting two different instanton solutions. In one case the 
probabilities are given by 



51x52 



,(a+VI^)^ 



(30) 



for 6 < ^4-, and in the other case the probabilities are 



Pg3 ~ e V ^ 



2A 



A2-f 1-5 



2A 



A2-f 1-5 



(31) 
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for < (5 < It is evident from eqs. (30) and (31) that the creation of a universe 

without PBH is more probable than with PBH for 5 < with a positive cosmological 
constant. It admits negative S also. However, the probability for creation of a universe 
with a pair of PBH is favoured in the first case when one considers a negative cosmological 
constant. In the second case, the probability of creation of a universe with a pair of PBH 
is favoured with a positive cosmological constant for < (5 < A^, but it is suppressed for 
A^ < (5 < ^ . In the second case no instanton solution exists with a negative cosmological 
constant. 

(d) For 6 = ^ and a ^ 0, the probabilities are given by 

Jr^gs ^ eA ' , J^giy^g2 ~ e3A ' . (^^2) 

In this case a positive cosmological constant and a > leads to a de Sitter universe creation 
more probable than a universe with a pair of black holes. However, for a < we note 
that creation of a universe with a PBH is more probable. 

(e) For a theory without cosmological constant, the evaluated probabilities are 

Fg3 ^ e V-3S ^ I-'gly^g2 ~ e 3V-35 _ 

It is evident that negative values for the parameter 6 are permitted to obtain a physically 
relevant solution. For 6 < 0, we note that a universe without a PBH is more probable for 
positive values of the prameter a, however, it is less probable if one considers a negative a 
satisfying the inequality a < ~i ^^3^ 1 determined by a negative S. 
IV. DISCUSSIONS : 

In this letter we estimated and compared the probability for creation of a universe with 
a primordial black hole (PBH) pair with a universe without such PBH pair considering a 
modified theory of gravity. In section II, we obtained gravitational instanton solutions in 
the two cases : (i) a universe with R x - topology and (ii) a universe with R x x 
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- topology respectively. The Euclidean action is then evaluated corresponding to the 
instanton solutions using BH technique [13]. We found that the probability of a universe 
with RxS^ topology turns out to be lower than a universe with topology RxS^xS^ in some 
specified restricted conditions. It may be mentioned here that one gets a regular instanton 
in 5^ - topology if there is no black hole. The existence of black hole restricts such a regular 
topology. The results obtained here on the probability of creation of a universe with a pair 
of primordial black holes are found to be strongly suppressed depending on the parameter 
S decided by the cosmological constant in some cases which are determined here. It is 
noted that the probability for a universe with PBH is strongly suppressed for small values 
of S lying in the sector < S < ^ in the presence of a positive cosmological constant. 
We note an interesting solution here in the framework of the modified gravitational action 
with inverse power of R -theory, which admits de Sitter instantons with and x 
topologies even without a cosmological constant. However, a negative 6 is required if one 
considers a theory without a cosmological constant. In this case the action is similar to 
that considered by Carroll et al. [6] to obtain an accelerating universe today. We also 
note another case where a de Sitter universe is less probable if a < — The result 
obtained earlier with curvature square term in the action in Ref. [17] given by eq. (29) 
now can be compared with eq. (32). It is evident that although the effect of inverse 
square term in the action enhances the probability, the condition to get a universe with a 
PBH remained unaltered in the theories here. It is interesting to note here that analytic 
continuation of a i? x 5^ metric considered here to Lorentzian region leads to an open 
3 - space. One obtains Hawking- Turok [23] type open inflationary universe in this case. 
In the other type of topology an open inflation section of the universe is not permitted. 
A detail study of an open inflationary universe will be discussed elsewhere. Thus in a 
modified Lagrangian with inverse power in i?-theory, quantum creation of PBH seems to 
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be suppressed in the minisuperspace cosmology for some values of the parameters in the 
action, which are determined here. Another new result obtained here is that gravitational 
instanton solution may be obtained even with a negative cosmological constant which is 
not permitted in the case considered by BH [13]. 
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